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I. Program Objectives 


We propose to perform spectroscopic studies, including data analysis and modeling, of star-formation 
regions using an ensemble of archival space-based data including some from the Infrared Space Ob- 
servatory’s Long Wavelength Spectrometer and Short Wavelength Spectrometer, and other spectro- 
scopic databases. We will include four kinds of regions: (1) disks around more evolved objects; (2) 
young, low or high mass pre-main sequence stars in star-formation regions; (3) star formation in ex- 
ternal, bright IR galaxies; and (4) the galactic center. During the proposed 2-1/2 year time frame of 
the program other data are likely to become available, for example perhaps from SIRTF, and we in- 
tend to be responsive to these and other such spectroscopic data sets. The OH lines in the far infra- 
red will be one key focus of this inquiry. 

The program has the following goals: 

1. Refine the data analysis of ISO observations, to obtain deeper and better SNR results on selected 
sources. The ISO data itself will be undergoing “pipeline 10” reductions in early 2001, and the 
more “hands-on data reduction packages” will be refined by the ISO teams into early/mid 2001. 
The Fabry-Perot database in particularly sensitive to noise and slight calibration errors, and im- 
provements are anticipated. We plan to build on these deep analysis tools, and contribute to 
their development; 

2. Model the atomic and molecular line shapes, in particular the OH lines, using revised Monte- 
Carlo techniques developed by the SWAS team at the Center for Astrophysics; 

3. Use newly acquired space-based or SOFIA spectroscopic data as they become available, and con- 
tribute to these observing programs as appropriate, as well as selected complementary radio 
and/or submillimeter results; 

4. Attend scientific meetings and workshops; 

5. Do E&PO activities related to infrared astrophysics and/or spectroscopy. 

II. Progress Report 

Work proceeded Fully on track and on time this past year. We review our accomplishments in the 
items below. 

A Workshops and Conferences 

One major item of the year program was participation in a planned and budgeted three-month work 
session at the Istituto Fisica Spaziale Interplanetario laboratory in Rome, Italy, where we worked 
with team collaborators on the FIR spectra obtained with the Infrared Space Observatory’s Long 
Wavelength Spectrometer. This Italian group is the headquarters for the ISO/LWS data on main 
sequence stars. This group is also active in the development of the Herschel spacecraft. One of the 
major outcomes of this effort are three conference papers, and one preprint, listed below (a copy of 
which is attached): 

“The Far Infrared Lines of OH as Molecular Cloud Diagnostics,” by Smith. H.A. . in the Proceedings 
of the 2 nd Maryland Conference on Far IR Astronomy, 2002. 

“IR Fine-Structure Line Signatures of Central Dust-Bounded Nebulae in Luminous Infrared Galax- 
ies,” by Fischer, J., Allen, R., Dudley, C., Satyapal, S., Luhman, M., Wolfire, W., and Smith, H.A-, in 
the Proceedings of the 2 nd Maryland Conference on Far IR Astronomy, 2002. 

“Mid and Far Infrared Spectroscopy of Seyfert Galaxies,” by Spinoglio, L., Malkan, M., Smith, H,A., 
Fischer, J., in the Proceedings of the Paris Conference on AGN, 2002 . 

“ISO/LWS Observations of the Bright Rimmed Globule IC 1396 N,” Saraceno, P., Smith, HtA», in 
preparation. 

A second major landmark was the hosting of a workshop at the CfA on “ISO/LWS Results on Ex- 
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tragalactic Astronomy from the Central Programme.” This meeting was held 6-7 December 2001 in 
Cambridge. It was very successful, and it appears likely that a follow-up workshop will be held this 
summer. 

Finally, as proposed, we spent time with our collaborator Dr. Glenn White, currently of the Univer- 
sity of Kent, Canterbury, England, working on the Galactic Center data set. We gave an invited 
seminar to the Physics Department there. 

B. Research Results 

During this period we successfully extended the SWAS (Submillimeter Wave Satellite) Monte-Carlo 
radiative transfer code to include: (1) OH; (2) the far IR lines of H 2 0; (3) [Oil; and are in processes of 
adding (4) the far infrared lines of CO. As a result we are able to model theoretically the emission 
and absorption from these species in molecular clouds whose structures we can adjust to ascertain 
the physical parameters most suitable for the various data observed. A fuller description is included 
in the attached preprint presented at the 2 nd Maryland Conference on Far IR Astronomy, 2002. 

Also during this period we began to use the “DUSTY” radiative transfer code to model the continuum 
emission from the same regions. The DUSTY code allows us to produce a fully thermodynamically 
self-consistent model of the cloud, from which we can then use the SWAS Monte-Carlo code to cal- 
culate the line emission properties. 


Also during this period we spent a considerable effort helping to prepare for SIRTF observations, 
including IRAC observations of extragalactic sources, and including the proposed Early Release Ob- 
servations (EROs). 

Finally, as expected, the ISO pipeline 10 data became available, and we have used these data as the 
basis for new reductions, and the analyses above. 

C. Advising Students and Postdocs 

During this past year the project PI has helped to advise one postdoctoral researcher, Dr. Lisa 
Kewley, who is now at the CfA. Her specialty is the AGN-IRB connection, and one conference pro- 
ceeding has so far been published: “Do Mergers Stop Monsters?” Kewley, L., Dopita, M., and Smith * 
HA.. . B.A.A.S., 199 , 4304, 2001. 

Also in this time period the PI advised the thesis research of two graduate students, who completed 
their theses: (1) Dr. Sarah Leeks, of Queen Mary and Westfield College, whose thesis was, “The 
Long Wavelength Spectrometer: Reduction and Interpretation of Data on W28 A2, a High-Mass 
Star-Forming Region”; and (2) Dr. Matt Bradford, of Cornell University, whose thesis was on Far 
Infrared Spectroscopy. 

D. Other Publications 

In addition to the publications listed above, we were co-author of the following: 

“An Infrared Study of the L1551 Star Formation Region - What We Have Learnt from ISO and the 
Promise for FIRST,” White, G. J., Liseau, Rene, Men’shchikov, A., Justtanont, K, Nisini, B., 
Benedettini, M., Caux, E., Correia, J. C., Kaufman, M., Molinari, S., Lorenzetti, D., Saraceno, P., 
Smith. H. A .. Spinoglio, L., Tommasi, E., and Fridlund, C. V. M., The Promise of the Herschel Space 
Observatory . Eds. G.L. Pilbratt, J. Cemicharo, A.M. Heras, T. Prusti, & R. Harris. ESA-SP 460, p. 
215, 2001. 

E. E&PO Activities 

During this period the PI gave several talks to elementary school groups. He also helped to prepare 
the E&PO program for NGST-based instruments. He spoke at the Boston Museum of Science to a 
assembly of museum staff on the upcoming cosmological “Cosmic Questions exhibition. 

III. Program Plans 

During the second year of this program we intend to continue apace with the program, as originally 
presented. We anticipate no particular difficulties. 
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IV. Reportable Patents/New Technology 


There have been no reportable patents/new technology made under this grant. 


Appendix A 

Preprint: “The Far Infrared Lines of OH as Molecular Cloud Diagnostics, by Smith, H,A; i to appear 
in the Proceedings of the 2 nd Maryland Conference on Far IR Astronomy, 2002. 
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To appear in the Proceedings of the 2 nd Workshop on New Concepts for Far-IR 

Submillimeter Space Astronomy, 2002 


- The Far Infrared Lines of OH as Molecular Cloud Diagnostics 

Howard A. Smith 

Harvard-Smithsonian Center for Astrophysics 
60 Garden Street 
Cambridge, MA 02138 
hsmith@cfa.harvard.edu 

Abstract 

Future IR missions should give some priority to high resolution spectroscopic 
observations of the set of far-IR transitions of OH. There are 15 far-IR lines 
arising between the lowest eight rotational levels of OH, and ISO detected 
nine of them. Furthermore, ISO found the OH lines, sometimes in emission 
and sometimes in absorption, in a wide variety of galactic and extragalactic 
objects ranging from AGB stars to molecular clouds to active galactic nuclei 
and ultra-luminous IR galaxies. The ISO/LWS Fabry-Perot resolved the 
1 19pm doublet line in a few of the strong sources. This set of OH lines 
provides a uniquely important diagnostic for many reasons: the lines span a 
wide wavelength range (28.9pm to 163.2pm); the transitions have fast 
radiative rates; the abundance of the species is relatively high; the IR 
continuum plays an important role as a pump; the contribution from shocks 
is relatively minor; and, not least, the powerful centimeter-wave radiation 
from OH allows comparison with radio and VLBI datasets. The problem is 
that the large number of sensitive free parameters, and the large optical 
depths of the strongest lines, make modeling the full set a difficult job. The 
SWAS montecarlo radiative transfer code has been used to analyze the 
ISO/LWS spectra of a number of objects with good success, including in both 
the lines and the FIR continuum; the DUSTY radiative transfer code was 
used to insure a self-consistent continuum. Other far IR lines including those 
from H 2 0, CO, and [OI] are also in the code. The OH lines all show features 
which future FIR spectrometers should be able to resolve, and which will 
enable further refinements in the details of each cloud’s structure. Some 
examples are given, including the case of S140, for which independent 
SWAS data found evidence for bulk flows. 

1 The Far Infrared Transitions of OH 

Storey, Watson and Townes (1981) made the first far infrared detection of OH in the 
interstellar medium: the two 1 19pm lambda-doubled lines between the ground and first 
excited states, which they discovered in absorption from Sag B2, and in emission from 
the shock in Orion KL. Altogether there are 15 far infrared lines between the lowest 8 
rotational levels of OH. The lines (each involving six transitions between the hyperfine- 
split rotational levels) are at wavelengths of (approximately) 28.9pm, 34.6pm, 43.9pm, 
48.8pm„ 53.3pm, 65.1pm, 71.2pm, 79.2pm, 84.4pm, 96.3pm, 98.7pm, 1 15.4pm, 
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1 19.9pm, 134.8pm, and 163.1pm. The upper state excitation temperatures for these 
lines range from 120K to 618K. The dipole moment for OH is large, 1 .668 Debye (for 
comparison the CO dipole moment is 0.1 12 Debye), and the radiative rates for OH 
transitions are generally fast. For example, the 1 19pm fundamental transition rate is 
about 0. 1 sec' 1 . But the FIR OH transitions also include some cross-ladder lines whose 
radiative rates are one hundred times weaker, providing a dataset of neighboring, far IR 
lines which frequently include both optically thin and very optically thick features. The 
OH analyses have an additional resource from which to draw: the strong hyperfine radio 
wavelength transitions that OH has in its ground-state, and which have been extensively 
observed. Maser and/or mega-maser activity is seen in many of the stars and galaxies 
which ISO observed. 

The general properties of interstellar OH are known from thermal OH emission studies 
done at radio wavelengths, as well as from the far IR observations. Typically OH has the 
following range of properties: N OH /N H2 = 0.1-3 xlO" 7 ; N OH L = 2-100 x 10 15 cm' 2 ; T x = 
100 - 275K; and N H2 = 0. 1 - 3 x 10 7 cm' 3 (e.g., Watson et al. 1985; Jones et al. 1994). In 
maser galaxies, such as the ones we observed with ISO/LWS, the OH masers regions 
typically have somewhat different properties: N 0H / N ra = 10' 7 - 10' 8 ; N OH L ~10 15 cm' 2 ; 
T x = 40K - 50K; N H2 ~ 0. 1-lxlO 4 cm' 3 (e.g., Henkel and Wilson, 1990). The presence of 
OH megamasers allow for VLBI observations in AGN, and results indicate that the sizes 
of the emitting regions are compact: ~ from a few to tens of parsecs. For example, in 
Arp220, one of the galaxies I discuss below, they seem to surround the AGN, with H 2 
densities ~ 10 6 cm' 3 (Lonsdale et al. 1994; Skinner et al. 1997). 

Four of the relatively strong far IR OH transitions involve the ground state. In warm 
clouds the molecules absorb the strong IR dust continuum, populating effectively some of 
the higher lying levels. As was conclusively shown by Sylvester et al. (1997), the 18 
centimeter radio maser emission in evolved stars is pumped by absorption of the 34pm 
dust continuum, while in AGN Skinner et al. (1997) proved the effectiveness of the IR 
pumping of mega-masers in Arp220. The presence of a strong far IR continuum affects 
all of the far IR lines, to varying degrees, and means (besides making the models more 
complicated) that the set of far IR OH lines also provide a sensitive measure of the local 
continuum conditions. One final point is worth noting about OH: unlike other some 
species commonly used as far IR probes of the interstellar medium, like CO, H 2 0, or OI, 
OH emission from PDRs or shocks does not contribute a relatively dominant amount of 
luminosity to these other processes. In the cases I model below, just the warm gas from a 
dusty molecular clouds is adequate to explain the strengths we observe. OH, therefore, 
has a powerful combination of features that makes it a very useful species for dis- 
entangling cloud properties across a very wide range of conditions. 
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2 ISO/LWS Observations of Extragalactic OH 

The ISO/LWS Extragalactic Science team has seen OH in fifteen galaxies, and has 
obtained potentially useful limits on approximately one hundred other galaxies observed 
by ISO. Our observations include a set on the 34.6pm “IR pump” transition from the OH 
2 PI 3/2 ground state, which was obtained using the ISO/SWS spectrometer. This line, and 
the less important 28.9pm line between excited upper states, are the only OH far IR lines 
not in the LWS wavelength coverage. ISO of course has also seen OH in numerous 
galactic locations, including Sag A, and Sag B2, and in particular in the evolved star 
IRC+ 10420 - the first source to demonstrate that the 34pm continuum can effectively 
pump the OH maser lines in these stars. The LWS extragalactic detections of OH are in 
the following sources: Arp 220, CenA, IRAS 17208-0014, IRAS201 00-4 156, M82, 
Mkn 231, Mkn 273, NGC 253, NGC891, NGC 1068, NGC 1614, NGC 3690A, 

NGC 4945, NGC 7469, and 3Zw35. Figure 1 shows all the eight OH lines detected in 
Arp220, along with a sample of the lines detected in other objects to illustrates some of 
the morphological variety. 

2.1 General Characteristics of the Observed Extragalactic OH Lines 
Perhaps the most striking observation about the set of extragalactic lines measured is the 
wide morphological range of behavior they display, even though all arise in infrared 
bright galaxies with either active star formation, an active nucleus, or perhaps both. In 
Arp 220, for example, every OH line is seen in absorption except the longest wavelength, 
163pm, line which is seen weakly in emission (Fischer et al. 1998) . By contrast, NGC 
1068, another AGN, has every detected line seen in emission, even the strong 1 19pm 
transition between ground and the first excited state (Spinoglio et al. 1999). NGC 253, a 
nearby starburst, has some OH lines in emission and some in absorption (Bradford et al. 
1999), while in M82, the infrared bright, prototype starburst galaxy, the lines’ equivalent 
widths are so small continuum that even with our high signal-to-noise ratio only the 
1 19pm line has been conclusively seen at all, in absorption (Colbert et al. 1999). From 
the analysis of this diverse set of line strengths, several useful preliminary generalizations 
may be drawn for the different categories, as summarized below. 

ULIGs: In the case of Arp 220, the VLBI megamaser studies together with the strength 
of the 34pm pump absorption provide some strong physical constraints. The OH lies in 
numerous small clouds which surround the AGN, and which are pumped by the local, 
warm, far IR continuum. Through further modeling I hope to arrive at a better sense of 
the shape of that continuum: is it starburst-like, or more AGN like, and is it characteristic 
of all ULIGs? 

AGN: In the case of the bright Seyfert galaxy, NGC 1068, analyses of the strong atomic 
lines (Spinoglio et al. 1999) show the substantial presence of a PDR line emission 
component, along with the high ionization lines from the AGN component in the ISO 
beam, and a starburst component, which is seen as well in other AGN. However the 
[CII]/FIR ratios are strange - often less than in PDRs. The OH lines, all seen in 
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emission, might help sort out the density and geometry of the clouds (for example, to see 
if the PDR regions have smaller than average filling factors) - and provide clues to their 
relationship to the active nucleus. 

Infrared Bright Galaxies: We find that the OH 1 19 pm ground fundamental transition is 
always in absorption in these galaxies, as is the 53pm line, while the 163pm feature is 
always seen in emission. Otherwise there does not appear to be any consistent behavior 
in the lines from different objects: in some sources they are seen in emission, in others 
they are in absorption. Saraceno et al. (1996) and Benedetti et al. (2000), among others, 
have noted there seems to be a dearth of H 2 0 emission in some galactic clouds. Also, 
the [CII] fluxes are very low, perhaps due to low gas heating (but this is not conclusively 
demonstrated). They note that the neutral oxygen [01] 63pm line is often abnormally 
weak and may be self absorbed. Finally, the SWAS satellite found that 0 2 is very weak 
or absent (Goldsmith et al. 2000). OH plays a key role in the chemistry of the ISM, is 
sensitive to the temperatures and radiation fields, and its abundance and distribution 
should help in the analyses of all these issues. 


3 Modeling of the Far Infrared OH Lines 

3,1 The Montecarlo Code for Lines: the “DUSTY” Code for Continuum 
I used the one-dimensional montecarlo radiative transfer code developed by the SWAS 
mission (Ashby et al. private comm) to model the OH line strengths. This code, a 
modification of the original Bemes code, adds a treatment of continuum photons from 
dust mixed in with the gas — a particularly essential feature for OH, which is pumped in 
many cases by absorption of 34pm continuum. In addition, the code corrects for some 
previous errors encountered at large optical depths, also an issue of importance for OH 
which has a very strong matrix element. Finally, the code includes an ability to handle a 
wide range of molecules besides rigid rotors. The montecarlo code takes as input a series 
of concentric shells, each of which is specified as to size, gas and dust temperature, H 2 
density, velocity and velocity width, and molecular abundance relative to H 2 . The model 
as a whole also assumes a (specifiable) dust emissivity. In all the modeling, the dust is 
assumed to be 1% of the gas, and to have the same temperature as the gas everywhere in 
the cloud. The montecarlo code calculates the populations of the molecular levels in 
each shell. The output is then fed into a radiative transfer code that calculates the line 
profiles as seen by an external observer looking at the cloud. 

The montecarlo results confirm the fact, noted above, that the OH lines are often optically 
thick. The peak optical depth at line center for the 1 19pm line exceeds 100 in some 
cases. As a result, many of the lines in many situations show self-absorption. As the 
model is tuned to fit the data by increasing the column density, these features can turn an 
emission line into an absorption line, one while strengthening the emission from a weaker 
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neighbor. The full set of nine line intensities, which span a range of wavelengths and 
optical depths, enable us to derive a rather detailed self-consistent picture of the cloud 
conditions solely from their intensity ratios, without the need for velocity information. 
But, as seen in the example below, there is also a wealth of information in the line shapes. 

ISO/LWS also observed lines of CO, H 2 0, and [OI] in these galaxies, with varying 
degrees of success. The SWAS montecarlo code can also predict the line emission from 
these species. In general the model predictions for these lines are in overall agreement 
with the observations; the uncertainties are due to uncertainties in the ISO line fluxes 
themselves, in the assumed molecular abundances, and in the amounts of possible line 
“contamination” from shocks and PDRs in the beam. Overall the results add further 
confidence to the models. The [OI] 63pm line is clearly self absorbed in several 
instances, confirming the suspicion that the line intensity is sometimes very weak due to 
self-absorption (Saraceno et al .) 

The montecarlo output is not strictly internally self-consistent; the input parameters need 
not conserve luminosity between shells, for example. To obtain this self consistency, and 
to guarantee that the final cloud model was also consistent with the observed far infrared 
continuum flux, I used the DUSTY code (Ivezic and Elitzur, 1997) to model the 
continuum and generate a set of shell parameters that provided this consistency. Then I 
iterated the DUSTY model with the montecarlo line outputs. While this technique does 
not give a unique solution for the cloud structure, it does provide a canonical model 
consistent with the observations. 


3.2 Some Modeling Results 

3.2.1 S 1 40 - A Molecular Cloud with Bulk Inflow Motions 

Ashby et al. (2000) used the SWAS satellite to observe S140, and the SWAS montecarlo 
code to model the shape of the observed submm H 2 0 line. They are able to obtain a 
detailed, if not entirely unique, model of the cloud. From their set of models they 
conclude, for example, that the cloud radius is 0.44pc, has an inner temperature of ~70K 
an inner hydrogen density of 1 .4x106 cm' 3 , a radial profile of temperature which varies 
like r° 5 , and a density profile varying like r' 08 . They also conclude that “significant bulk 
flow” is required to explain the H 2 0 line shape, but because of the small optical depth of 
the 557 GHz line they could not differentiate between infall and outflow. Although the 
ISO observations of SI 40 were only able to set weak limits on the OH lines (Aannestad 
and Emery, 1998;, the detailed nature of the SWAS models made it a useful check. 

Taking the SWAS cloud parameters for the case of bulk infall motions, I used the 
montecarlo code to calculate the strengths and line shapes for the full set of OH lines; I 
calculated the far IR lines of H 2 0 and the [OI] lines as well, for comparison with the ISO 
values. The reasonable results obtained from the model of SI 40, as confirmed with the 
SWAS (and ISO) observations, provide some confidence in the galaxy modeling. Figure 
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2 shows the results for three of the OH lines, on the same scale: the 163pm line (top), the 
79pm line (middle), and the 1 19pm line (bottom). It is clear that high spectral resolution 
observations of the far IR water lines can readily distinguish infall from outflow, because 
these lines are optically thicker than the submillimeter line. 

3.2.2 Arp220 - The Molecular Cloud Component of a Peculiar Ult ra Luminous Galaxy. 
Arp220 is a particular challenge, because so many OH lines are seen, and every one of 
them (except the line at 163pm) is seen in absorption — something that happens in no 
other known object. Arp220 is also unusual in general in that its spectrum is charac- 
terized by numerous molecular absorptions; even the [OI] 63pm fine structure line is seen 
in absorption, while the [CII] 157pm line is found only in weak emission (e.g., Fischer et 
al. 1998). Skinner et al. (1997) showed that the 34pm OH absorption could pump the 
OH megamasers in Arp220, and their analysis supported the model in which many small, 
dense molecular clouds circulate around the nucleus. Suter et al. (1998) attempted to 
model the OH line absorptions with some less complex radiative transfer models, and 
were driven to consider some unusual non-LTE situations to explain the observations; . 

I used the SWAS montecarlo code successfully to model (to first order) all of the 
observed OH lines in Arp220. Figure 3 shows a sample of these model line shapes. In 
order to get absorption in all the lines, including the 79pm line to the ground state PI3/2 
from the lowest level of the 2 PI 1/2 ladder, while still getting weak emisson at 163 pm, it 
was necessary to have falling within the ISO/LWS beam a combination small dense 
clouds, and a few large, giant molecular clouds. Iterating with the DUSTY code provides 
a way to obtain a continuum shape that fits the observed infrared continuum, though there 
may also be a continuum component in our beam arising from clouds without much OH. 
It was not necessary in this process to resort to the non-LTE scenario postulated in Suter 
et al. 

3.3.3 NGC 253 - The OH Emission from a Nearby Infrared Bright Galaxy 

In NGC 253, the two strongest OH lines are seen as emission lines: the 79pm and 163pm 
features, each of which is about twice as strong as the 1 19pm fundamental absorption 
(Bradford et al. 1999). Only two other OH lines are also seen in this galaxy, and as a 
result the model’s constraints are weaker than for Arp220. The montecarlo modeling of 
these OH lines implies that a few giant molecular clouds, with r ~ lOOpc, can explain the 
observations. The clouds have H 2 densities averaging approximately 3x1 0 3 cm 3 and an 
OH abundance relative to H of about 5x1 O' 8 . ISO/LWS also observed the 1 19pm line in 
NGC 253 with the high resolution Fabry-Perot, which Bradford (2001) has successfully 
modeled. 

3.3.4 NGC 1068 - OH Emission from AGN 

NGC 1068 is a remarkable extragalactic source in its OH spectrum, because it is the only 
galaxy observed in which all the detected transitions are found in emission. The code 
accurately predicts the observed line flux ratios to 40%. Based on the montecarlo 
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modeling, it seems the clouds that produce these lines in NGC1068 must be relatively 
small and dense, and heated from the inside. They are about 0.2pc in radius, with 
densities of < 1 0 4 cm' 3 at the outer edges, increasing towards the center with a power law 
behavior of (R/Ro ) 1 25 . The temperature in the outer shells is about 25K, increasing 
inwards approximately with a power law dependence (R/Ro ) 1 41 . The total column 
density of H 2 in each cloud is 1 .5xl0 24 cm' 2 , and the relative abundance of OH in the 
model is approximately 1 O' 7 . OH is taken to be absent in the hotter portions of the cloud, 
where T>300K. With these high column densities the strong lines of OH are very 
optically thick, and precise radiative transfer calculations like the present ones are 
absolutely essential. At a distance of 1 6.2Mpc, the total number of equivalent number of 
clouds needed to produce the observed absolute flux in the lines and the continuum is 
~3xl0 7 . The limits to the observed fluxes in all the weaker lines are consistent with this 
modeling. It is worth noting in this context that in some galaxies the weaker OH lines in 
NGC 1068, like the 53pm line, are amongst the strongest, but the montecarlo modeling 
can successfully account for these differences. The montecarlo code also successfully 
reproduces the continuum emission seen from NGC 1068 to within a factor of 2-3 in 
absolute flux density across the entire LWS spectrum. The total mass in such an ensemble 
of clouds in NGC 1068 is approximately 9xl0 9 Mq. 


4 Conclusions to Date 

What general conclusions might we hope to draw from this large set of ISO OH 
observations? Although the analysis is still underway, it appears that the infrared 
luminous galaxies can be grouped into three general categories based on the relative 
strengths of their infrared OH lines: (1) relatively normal galaxies like NGC 253 and 
M82, from which the OH lines are seen in both emission and absorption, and which are 
dominated by starburst activity in giant molecular clouds; (2) NGC 1068 and other AGN, 
whose OH is seen in emission, and whose clouds must be quite small and dense 
(hydrogen densities up to 10 s cm" 3 in the cloud cores); and (3) Arp220, and perhaps other 
ULIGs with active starbursts as well as an active nucleus, for which a combination of a 
few giant molecular clouds and a collection of moderately dense cloudlets are required. 
(The ISO/LWS spectrum of Mkn 231 is much noisier than Arp 220’s, but hints at similar 
behavior in its OH.) While I have made substantial progress using only the intensities of 
these lines, velocity resolved spectra will enable even further refinements. Future far 
infrared missions with sensitive, spectroscopic capabilities should find that the set of OH 
lines can provide wealth of information needed to unravel the structures of complex 
clouds. 
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Appendix A 



Sample of next four lines: 



Figure 1 : The 8 ISO/LWS OH lines seen in Arp 220 (the 34pm SWS-band line is not 
shown here). Under each OH line in Arp220, an OH line selected from another other 
ISO galaxy (or IRC+ 10420) is shown for comparison at the corresponding wavelength to 
illustrate some of the variety in the intensities observed. 
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Figure 2: The hydrogen density distribution in a nominal giant molecular cloud in the 
galaxy NGC 253. This distribution was used in the montecarlo code modeling. The 
inner 50 parsecs of the GMC has warm, 250K, gas and dust, at constant density, while the 
bulk of the cloud has much cooler, 35K material. 



Figure 3: Predicted OH line shapes from infalling gas in SI 40, based on cloud profiles 
derived from SWAS data. The figure shows the 163 pm line (top), the 79 pm line 
(middle), and the 1 19 pm line (bottom). While all three lines show evidence for infall 
from their highly asymmetric shapes, only the latter two show characteristic absorption: 
in the redshifted (near side) material. The velocity range across the figure is ±10 km/sec, 
reflecting the fact that the maximum infall velocity in the model (and seen by SWAS) 
was 7.2 km/sec. 
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Velocity 



Figure 4: Montecarlo modeling results of three OH lines in Arp 220. The OH lines at 
1 19pm, 79 pm, and 53 pm are displayed. All the OH lines are seen in absorption in Arp 
220, except the weakly emitting 163 pm line - the only known object with this property. 
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